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Preface 


This  thesis  is  a  continuation  of  work  previously  done  at  the  Air 
Force  Institute  of  Technology  using  Monte  Carlo  sampling  techniques  to 
conduct  sequential  probability  ratio  tests  of  the  Weibull  density 
function.  It  is  hoped  that  this  thesis  will  provide  more  accurate  and 
dependable  test  plans  of  reliability  when  the  underlying  distribution 
is  a  two-parameter  Weibull. 

I  would  like  to  thank  ray  advisor.  Dr.  Albert  H.  Moore,  for  both 
suggesting  ray  topic,  and  providing  guidance  as  the  work  progressed. 


ii 


Table  of  Contents 


Preface  . 

List  of  Figures  . 

List  of  Tables  . 

Abstract  . 

I.  Introduction  . 

Purpose  . 

Verification  . 

II.  Background  . 

The  Wei  bull  Distribution . 

Sequential  Tests  of  Hypotheses  . 

Sequential  Tests  with  the  Weibull  Distribution  . 

III.  Methodology  . 

Assumptions  . 

Generation  of  Random  Weibull  Deviates  . 

Maximum  Likelihood  Estimation  of  Weibull  Parameters  .  .  .  . 
Determination  of  Variances  of  R  by  Asymptotic 

and  Empirical  Means  . 

Sequential  Probability  Ratio  Tests  of  Reliability  .  .  .  . 
Computerization  . 

IV.  Results  and  Recommendations  . 

Results  . 

Recommendations  . 

Bibliography  . 

Appendix  A:  Computer  Programs  . 

Appendix  B:  Tables  . 

Vita . 


iii 


List  of  FI qures 


Figure 

1.  Weibull  Density  Function  for  9=i  .  .  . 

2.  Graphical  Representation  of  a  SPRT  .  . 

3.  Truncated  Log  SPRT  . 


iv 


mmm 


List  of  Tables 


Table  Page 

B-1  Sequential  Test  Results . 41 

B-2  Sensitivity  Analysis  .  42 

B-3  Sequential  Test  Results  with  Shift . 43 

B-4  Sequential  Test  Plans  with  Shift  and 

Revised  Input  Risk  Levels  .  44 


V 


AFIT/G0R/MA/81D-8 


Abstract 


Monte  Carlo  Analysis  techniques  are  used  for  the  development  of 
test  plans  concerning  hypothesized  system  reliabilities.  Systems  under 
consideration  are  those  in  which  component  failure  rates  are  best 
described  by  the  two-parameter  Weibull  probability  density  function. 

The  statistical  test  employed  is  Wald's  sequential  probability  ratio 
test  using  independent,  asymptotically  computed  variances  (Cramer-Rao 
Lower  Bound  Technique),  The  null  hypothesis,  HqI  Rq  =  .90,  is  tested 
against  alternative  reliabilities  of  .854,  .810,  and  .729.  Three  pairs 
of  alpha  and  beta  risk  levels  are  considered  for  each  test  ((.l,.l), 
(.2,. 2),  (,1,,2)),  A  truncation  decision  for  the  sequential  test  is 
made  at  1.5  times  the  fixed  sample  test  size  for  the  same  conditions. 

One  thousand  Monte  Carlo  repetitions  are  used  for  these  test  procedures. 


Sequential  Testing  of  Hypotheses  Concerning 
the  Reliability  of  a  System  Modeled  by  a 
Two-Parameter  Wei  bull  Distribution 


I .  Introduction 

The  importance  of  reliability  testing  is  reflected  in  the  constantly 
increasing  emphasis  placed  on  this  subject  by  both  government  and  com¬ 
mercial  industry.  Most  Department  of  Defense  contracts  impose  some 
degree  of  reliability  requirements  on  the  contractor.  These  range  from 
the  definition  of  system  reliability  goals  to  the  requirement  for  actual 
demonstration  of  achievement.  Many  of  these  have  specific  funds  allo¬ 
cated  to  the  reliability  effort.  Some  also  require  the  development  and 
maintenance  of  a  reliability  program  plan,  and  specify  the  preparation 
of  periodic  reliability  reports.  Statistical  testing,  and  in  particular 
Monte  Carlo  sampling  techniques,  have  proven  quite  useful  when  testing 
and  evaluating  the  reliability  of  a  system. 

Purpose 

The  purpose  of  this  thesis  is  to  develop  truncated  (truncation 
occurs  at  1.5  times  the  fixed  sample  test  size)  sequential  test  plans, 
concerning  hypothesized  system  reliabilities.  Testing  applies  to  those 
systems  that  can  be  modeled  by  the  two-parameter  Wei  bull  probability 
distribution.  The  test  used  is  Wald's  Sequential  Probability  Ratio 
Test  (SPRT).  Wald  states  that  the  SPRT  "frequently  results  in  a  savings 
of  about  50  percent  in  the  number  of  observations  over  the  most  efficient 


test  procedure  based  on  a  fixed  number  of  observations"  (Ref  28:1). 

The  null  hypothesis,  Hg:R  =  .90,  will  be  tested  against  alternative 
reliabilities  of  .854,  .810,  and  .729.  Three  sets  of  a  and  B  risk 
levels  will  be  used  as  input:  (.2, .2),  (.l,.l),  and  (.1,.2). 

Verification 

Actual  alpha  errors  will  be  evaluated  by  inputting  Hq  true,  while 
actual  beta  errors  can  be  determined  by  inputting  true.  Verification 
of  these  test  plans  will  include  evaluation  of  these  Type  I  and  Type  II 
errors  using  Monte  Carlo  simulation.  Previous  work  at  the  Air  Force 
Institute  of  Technology  (Ackerson  1977,  Ballard  1978,  Antoon  1979)  have 
shown  Monte  Carlo  techniques  quite  useful  in  Wei  bull  testing. 

The  test  statistic  employed  for  the  sequential  tests  will  be  a 
modification  of  the  one  used  by  Ballard  in  his  1978  thesis  work.  The 
statistic  is  based  on  a  Tiklihood  ratio  of  the  hypothesized  reliabili¬ 
ties,  Rg  and  R^. 


II.  Background 


The  Wei  bull  Distribution 

In  the  late  1930' s,  the  subject  of  fatigue  life  in  materials  and 
the  related  subject  of  extreme  value  theory  were  being  studied  by 
Weibull  (Ref  30),  Gumbel  (Ref  8),  and  Epstein  (Ref  6)  among  others.  In 
1939  Waloddi  Weibull,  a  professor  at  the  Royal  Institute  of  Technology 
in  Sweden,  proposed  the  distribution  named  after  him  as  an  appropriate 
distribution  to  describe  the  life  length  of  materials.  In  1951,  while 
working  for  the  A.B.  Bofors  Steel  Company  in  Stockholm,  Sweden,  Weibull 
demonstrated  that  this  function  could  be  used  to  model  such  things  as; 
yield  strength  of  Bofors  steel,  size  distribution  of  fly  ash,  fiber 
strength  of  Indian  cotton,  length  of  Cytroidae,  fatigue  life  of  steel, 
statures  for  adult  males  born  in  the  British  Isles,  and  breadth  of  beans 
of  Phaseolus  Vulgaris  (Ref  29). 

Of  interest  to  the  Air  Force,  is  the  fact  that  the  Weibull  density 
function  can  be  used  to  model  aircraft  subsystem  life,  the  failure  and 
fatigue  life  of  metals,  the  life  of  many  electrical  components,  as  well 
as  many  other  items  currently  in  Air  Force  inventories.  The  probability 
density  function  for  the  three-parameter  Weibull  is  given  by: 

f(x;9,k,c)  =  k(x-c)*^~^  exp  C-((x-c)/9)*^]  x^o  (1) 

9*< 

®  9,k>o 

=  0  elsewhere 
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The  three  parameters  are; 
k-the  shape  parameter 
9-the  scale  parameter 
c-the  location  parameter 

The  shape  parameter,  k,  determines  the  shape  or  amount  of  peak  to 
the  curve.  This  parameter  allows  the  failure  rates  to  decrease,  in¬ 
crease,  or  remain  constant  over  time.  In  reliability  testing,  k  is 
considered  a  function  of  mean  ultimate  strength  (Ref  9:32)  (note  that 
for  k  equal  to  one  the  distribution  becomes  exponential). 

The  scale  parameter,  e,  is  sometimes  referred  to  as  the  charac¬ 
teristic  life  of  the  component  being  tested,  and  determines  the  spread 
of  the  function  about  it's  mean.  In  reliability  testing,  e  is  con¬ 
sidered  a  function  of  stress  (Ref  9:32). 

The  location  parameter,  c,  is  the  value  of  x  at  which  failures 
begin  to  occur.  In  this  thesis  c  will  be  zero,  reducing  the  three- 
parameter  Weibull  distribution  to  two  parameters.  This  indicates  that 
failures  can  occur  immediately  after  initiation  of  an  operation  or 
function.  The  two-parameter  Weibull  distribution  has  the  following 
characteristics: 

1.  Probability  density  function  (p.d.f.): 

f(x;0,k)  =  k(x)*^~^  exp  [-(x/e)*^]  x>o  (2) 

0*^  . 

9  ,k>o 

=  0  elsewhere 
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2.  Cumulative  distribution  function  (c.d.f.): 

F(x)  =  1  -  exp  [-(x/e)*^]  x>o  (3) 

e,k>o 

=  0  elsewhere 

3.  Reliability  function  (1  -  F(x)): 

R(x)  =  exp  [-(x/9)*^]  x^o  (4) 

8  ,k>o 

=  0  elsewhere 


The  range  of  shapes  that  a  graph  of  the  Wei  bull  density  function 
can  take  on  is  very  broad,  depending  on  the  value  of  the  scale  para¬ 
meter  9,  and  the  shape  parameter  k.  Figure  1  (Ref  4:31)  shows  three 
of  those  curves  with  9=1,  corresponding  to  k=l/2,  k=l ,  and  k=3. 


Sequential  Tests  of  Hypotheses 

Sequential  tests  of  hypotheses  differ  from  fixed  tests  in  that  the 
sample  size  is  a  random  variable.  After  each  trial,  one  of  three 


decisions  must  be  made: 


1.  Accept  the  null  hypothesis 

2.  Reject  the  null  hypothesis 

3.  Take  another  observation  and  continue  testing 

The  decision  made  is  based  upon  a  probability  ratio  test  statistic. 
As  mentioned  previously  this  testing  procedure  was  first  developed 
by  Wald. 
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Figure  1.  Weibull  Density  Function  for  9=1 


In  any  statistical  test  of  hypothesis,  there  exists  three  possible 
results:  the  correct  decision,  a  Type  I  error,  or  a  Type  II  error. 

The  correct  decision  is  when  the  test  fails  to  reject  the  null  hypo¬ 
thesis,  when  it  is  true.  A  Type  I  error  is  when  the  test  rejects  the 
null  hypothesis,  when  it  is  true.  A  Type  II  error  occurs  when  the  test 
fails  to  reject  the  null  hypothesis,  when  the  alternative  hypothesis  is 
true.  These  two  errors  are  referred  to  as  the  Alpha  and  Beta  errors, 
where  Alpha  is  the  probability  of  a  Type  I  error  occurring,  and  Beta 


is  the  probability  of  a  Type  II  error  occurring.  Alpha  is  sometimes 
called  the  producer's  risk,  while  Beta  is  often  referred  to  as  the  con¬ 
sumer's  risk  (Ref  7:48). 

Wald's  sequential  probability  ratio  test  considers  only  two  cases: 
1)  X  admits  a  probability  density  function;  2)  x  has  a  discrete  dis¬ 
tribution  (Ref  28:37).  In  other  words,  for  the  continuous  case  (such 
as  the  Weibull)  the  density  function  must  be  known.  The  three  decision 
regions  (acceptance,  rejection,  continue  testing)  were  successfully 
defined  by  Wald,  provided  the  test  statistic  is  based  upon  a  probability 
ratio.  Upper  and  lower  bounds  are  constructed  using  the  desired  Alpha 
and  Beta  error  levels  (Ref  28:40-42). 

1 .  Upper  Bound 

A  =  (1  -  Beta)/Alpha  (5) 

2.  Lower  Bound 

B  =  Beta/(1  -  Alpha)  (6) 

The  decision  rules  are: 

1 .  Accept  Hq  if  £  B 

2.  Reject  Hq  if  £  A 

3.  Take  another  observation  and  continue  testing 

if  B  <  <  A 

(note:  Z^  is  equal  to  the  likelihood  ratio  or  test  statistic) 
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Figure  2  is  a  graphical  representation  of  the  sequential  probability 
ratio  test  (Ref  7:85),  in  which  the  decision  to  accept  Hq  is  made  on  the 
eighth  observation. 


Rejection  region 


* 

* 

■k 

★ 

1 

■  CVJ 

3 

4 

★ 

5 

* 

6 

7 

* 

00  . 

Acceptance  region 


Figure  2.  Graphical  Representation  of  a  SPRT 


Sequential  Tests  with  the  Wei  bull  Distribution 

Work  on  SPRTs  with  the  Wei  bull  distribution  began  in  the  early 
1970s.  In  1971,  Nicolae  and  Obreja  developed  a  sequential  test  for  the 
two- parameter  Wei  bull  distribution  with  known  shape  parameter  (Ref 
20:320-331).  Since  that  time,  additional  work  has  been  done  at  the 
Air  Force  Institute  of  Technology.  In  1975,  Callahan  derived  formulas 
for  determining  values  for  the  test  statistic,  the  limits  A  and  B,  the 
expected  time  to  failure,  and  the  expected  number  of  failures  to  a 
decision.  His  work  concerned  both  discrete  and  continuous  samples 
for  SPRTs,  accelerated  SPRTs,  and  truncated  SPRTs  for  the  two-parameter 
Wei  bull  distribution  with  known  shape  parameter  (Ref  5).  Williams 
conducted  SPRTs  on  the  two-parameter  Wei  bull,  with  known  shape  para- 
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meter,  to  determine  the  effects  of  truncation  points- (Ref  31). 

Robinson  developed  a  standardized  set  of  SPRTs  for  use  with  the  Wei  bull 
distribution  when  the  shape  parameter  is  known  (Ref  24).  Hoffert 
studied  composite  SPRTs  for  the  Weibull  scale  parameter  with  the  shape 
parameter  unknown  (Ref  13).  Ackerson  (Ref  1)  and  Ballard  (Ref  3)  de¬ 
veloped  new  SPRTs  for  the  two-parameter  Weibull  distribution  using 
maximum  likelihood  estimates  for  the  scale  and  shape  parameters. 

Antoon  developed  different  methods  to  compute  the  variance  of  different 
reliability  estimates  for  systems  modeled  by  the  two-parameter  Weibull 
(Ref  2).  Monte  Carlo  Analysis  Techniques  were  used  extensively  in 
these  theses. 


Ill .  Methodolo' 


The  objective  of  statistics  is  to  make  inferences  about  a 
population  based  on  the  information  contained  within  the  sample  data 
(Ref  17:325).  Inferences  can  be  made  by  either  testing  hypotheses  about 
population  parameters  or  by  estimating  the  population  parameters.  A 
statistical  hypothesis  is  a  statement  about  the  parameters  of  the  random 
variable{s).  A  simple  hypothesis  completely  describes  the  distribution 
while  a  composite  hypothesis  does  not.  For  example,  a  hypothesis  of  the 
form  H:  a  £  ID  is  called  simple  if  (d  consists  of  a  single  point,  and 
composite  if  id  has  more  than  one  element  (Ref  18:281  ,290). 

The  reliability  of  a  system  is  defined  to  be  the  probability  that 
the  system  will  still  be  operating  under  stated  environmental  conditions 
at  some  specified  time  t  (Ref  16:8).  This  thesis  will  be  concerned 
with  this  reliability  based  on  failure  rates  that  start  from  initial 
component  use  until  failure  of  the  component. 

■Assumptions 

1.  Both  the  null  and  alternative  hypotheses  are  simple,  and  of 


2.  Failure  rates  are  best  described  by  the  two-pararaeter  Weibull 
probability  density  function  which  yields  the  following  reliability 
function: 

R{t)  =  exp  (-(t/e)*^)  t^o  C7) 

k,0>o 

=  0  elsewhere 

where  R  is  the  reliability,  t  the  time,  and  k  and  e  are  the  shape  and 
scale  parameters  of  the  Weibull  distribution  (Ref  25:72), 

3.  The  distribution  of  the  estimated  reliabilities  is  assumed  to 
be  Normal  and  dependant  only  upon  the  true  reliability,  R,  and  the  sampl 
size,  n. 

Generation  of  Random  Weibull  Deviates 

Random  Weibull  deviates  were  generated  using  the  CDC  6600  computer 
and  the  International  Mathematical  and  Statistics  Library  (IMSL)  sub¬ 
routine  GGEXN  (Ref  14),  to  first  generate  exponential  deviates.  The 
mean  and  standard  deviation  is  required  as  input  to  GGEXN.  Recall  that 
for  k  equal  to  1,  the  Weibull  distribution  becomes  an  exponential  dis¬ 
tribution  with  mean  and  standard  deviation  both  equal  to  9.  Solving 
equation  (7)  for  9  yields: 

=  -(t'^)/ln  R(t)  (8) 

t( 

This  value  of  9  was  used  as  input  to  GGEXN  to  create  random  deviates 
with  the  hypothesized  reliabilities.  These  deviates  are  then  raised  to 
the  1/k  power  which  yields  a  random  Weibull  deviate  (Ref  12:406).  The 
value,  R(t),  is  referred  to  as  the  input  stream  (Ref  3:10).  When 
generating  Weibull  deviates,  if  Rq  is  used  for  R(t)  the  data  generated 
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is  used  to  test  the  null  hypothesis  (input  Hq  true  to  determine  the 
probability  of  a  Type  I  error.  Alpha).  When  is  used  for  R(t)  the 
data  generated  is  used  to  test  the  alternate  hypothesis  (input  true 
to  determine  the  probability  of  a  Type  II  error,  Beta). 

Maximum  Likelihood  Estimation  of  Wei  bull  Parameters 

The  procedure  used  to  derive  the  maximum  likelihood  estimates 

A 

(MLE),  k  and  0,  of  the  Wei  bull  shape  and  scale  parameters,  k  and  0, 
was  developed  by  Harter  and  Moore  (Ref  11).  They  developed  a  procedure 
for  censored  or  uncensored  samples  for  the  three-parameter  Wei  bull  dis¬ 
tribution.  The  likelihood  function  for  the  Weibull  distribution  is: 

n 

LCx,c,0,k)  =  n  pdf  (x^:c,0,k)  (9) 

i=l 

The  natural  logarithm  of  the  likelihood  function  is: 

In  L  *  In  n!  -  In  (n-m)!  -  In  r!  +  (m-r)(ln  k-k  In  e) 

+  (k-1)  i  In  (x^-c)  -  L  [(x^-c)/0]^ 
i=r+l  i=r+l 

-  Cn-m)  [Cx^-c)/eJ*^  +  r  In  {1  -  exp  [-(x^.-c)/©]*^}  (10) 

The  first  partial  derivative  of  equation  (10),  with  respect  to  0  is: 

k  k+1 

L-  =  -k(m-r)/9  +  k  Z  (x.-c)  /0^  +  k(n-m) 

9  1 

i=r+l 

(Xi^-c)*^/©*^"^^  -  kr  (x^^^-c)*'  exp  [-(x^^^-c)'^/©'^] 

{1  -  exp  [-(x^^^^c)"/©'']}  (11) 
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The  first  partial  derivative  of  equation  (10),  with  respect  to  k  is: 

l-|^  =  (m-r)(1/k-ln  0)  +  e  In  (x^-c)  -  E  [(x^.-c)/0]^ 

i=r+l  i=r+l 

In  [(x^.-c)/0]  -  (n-m)  [(x^-c)/0]*^  In  [(x^-c)/0] 

+  ®^P  {-[(x^+Tc)/0]'^} 

70''  {1  -  exp  C-(x^^^-c)''/0'']}  (12) 

where  the  x^  are  the  Weibull  deviates,  n  is  sample  size,  m  is  the  first 
order  statistics  of  sample  size  n,  r  is  the  number  of  deviates  censored 
from  below,  and  (n-m)  is  the  number  of  deviates  censored  from  above 
(Ref  11).  Since  this  thesis  deals  with  the  two-parameter  Weibull  dis¬ 
tribution,  as  opposed  to  three,  and  there  is  no  censoring,  c=0,  (n-m)=0, 
and  r®0.  Substituting  these  values  into  equations  (10),  (11),  and  (12) 
yeilds: 

n 

In  L  =  In  n!  +  n  (In  k-k  In  0)  +  (k-1)  e  In  x^. 

i=l 

"  k 

-  I  (x./9)  (13) 

i=l  ^ 

Lg  =  -  k  n/0  +  k  s  (x.'^/0'^'^^)  (14) 

i=l 

n  n  . 

L|^  =  n/k  -  n  In  0  +  E  In  x^  -  E  (x^/0)  In  (x^./0)  (15) 

i=l  i=l 
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Equations  (.14)  and  (15)  are  then  set  equal  to  zero  and  solved  for  9 
and  k  respectively.  Equation  (14)  set  equal  to  zero  can  be  solved 
directly  for  9  given  k: 

e  =  (I  x.^/n)^/''  (16) 

i*l 

This  expression  is  then  substituted  into  equation  (15)  for  6.  Now, 
since  equation  (15)  set  equal  to  zero  is  simply  a  function  of  the 
deviates,  x..'s,  sample  size,  n,  and  k,  a  root  to  this  equation  can  be 
approximated  iteratively,  using  any  one  of  a  number  of  different 
numerical  methods.  The  method  used  here  is  the  bisection  method  (Ref  14), 
and  9  is  reevaluated  at  each  successive  value  of  k, 

Thoman,  Baine,  and  Antle  have  provided  an  unbiasing  factor  for  the 
MLE  of  k  (Ref  26).  However,  Petrich  showed  that  as  sample  size  increases, 
this  factor  approaches  one  and  can  be  neglected  for  samples  as  small  as 
six  (Ref  22:28-30).  The  unbiasing  factor  was  not  used  in  this  thesis 
since  the  minimum  sample  size  is  ten. 

Other  possible  methods  of  estimating  k  and  9  are  the  linear 
regression  method  (Ref  21:48-50)  and  the  method  of  matching  moments 
(Ref  23). 

Once  k  and  9  have  been  determined,  an  estimate  for  the  reliability 
may  be  found  using  the  formula: 

R  =  exp  [-(t/e)*^]  t^o  (17) 

k  ,9>o 

The  estimate  for  the  reliability,  R,  is  essential  for  computing  the 
variance  and  the  subsequent  sequential  tests. 


I  ‘i  I  **  ■ifii 
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Detennination  of  Variance  of  R  by  Asymptotic  and  Empirical  Means 

Empirical  methods  for  determining  the  variance  of  point  estimators 
are  given  by  Mendenhall  and  Scheaffer  (Ref  17:269),  for  both  biased 
and  unbiased  estimations.  In  addition  to  this,  Antoon  developed  a  pro¬ 
cedure  for  empirically  computing  the  variance  of  R  using  a  program 
called  FITIT  (Ref  2:22-25).  Ballard,  in  his  sequential  tests,  uses  a 
cubic  approximation  to  empirically  determine  the  variance  of  R.  He 
also  computed  variances  asymptotically  using  equations  from  Air  Force 
technical  training  notes  (Ref  3:13-14). 

In  this  thesis  the  variance  of  R  used  in  the  sequential  testing 
is  computed  by  the  Cramer-Rao  Lower  Bound  (CRLB).  Thoman,  Bain,  and 
Antle  (Ref  27)  determined  that  R,  the  maximum  likelihood  estimator  of 
R,  the  true  reliability  (assuming  a  two-parameter  Weibull  distribution), 
is  very  nearly  unbiased  and  has  a  variance  that  is  very  nearly  equal  to 
the  CRLB  for  the  variance  of  an  unbiased  estimator.  The  CRLB  is  a 
function  of  R  and  n,  and  is  given  by  the  following  formula: 

CRLB  =  R^/n  (In  R)^  [1.109  -  .514  In  (-In  R)  +  .608  (In  (-  In  R))^] 

(18) 

Here  it  is  assumed  that  the  CRLB  is  an  asymptotic  method  for  calculating 
the  variance  of  the  estimated  reliabilities. 

Sequential  Probability  Ratio  Tests  of  Reliability 

Previous  SPRTs  of  reliability  at  the  Air  Force  Institute  of 
Technology  have  been  conducted  by  Ackerson,  Ballard,  and  Jewell.  This 
thesis  is  an  attempt  to  extract  the  most  successful  methodologies, 
along  with  the  most  recent  test  statistic,  to  provide  more  dependable 
and  accurate  test  plans. 
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Ackerson  formed  the  following  test  statistic  based  on  the 


likelihood  ratio  of 


=  In  (Rq/R^)  +  R^  (-R^  +  2^  -  R^^) 

o  p  2  2^ 

2  R^  0  p 


+  R^  (R^  -  2RRq  +  Rq^) 

K »  'r 


(19) 


where  R  is  the  reliability  as  estimated  from  the  samples  using  MLE  of 
k  and  0,  found  by  placing  k  and  0  in  equation  (7).  Rq  is  the  relia¬ 
bility  under  the  null  hypothesis  and  R^  is  the  reliability  under  the 
alternate  hypothesis.  Ackerson  computed  the  variance  of  the  estimated 
reliabilities  using  sets  of  equations  from  the  Air  Force  technical 
training  notes. 

Ballard  developed  and  used  the  following  test  statistic  for  his 
sequential  tests: 


\  =  In  (L(R^)/L(Rq))  =  2RR^  -  2RRq  +  Rq^ 


2  0 


2. 

R 


(20) 


where  the  parameters  are  defined  the  same  with  the  exception  of  the 
variance  calculation,  Ballard  computed  the  variance  of  the  estimated 
reliabilities  using  a  cubic  approximation  generated  by  empirical 
means. 
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The  test  statistic  employed  in  these  sequential  tests  is  formed 


from  the  following  likelihood  ratio: 


L(R■^) 

KRq) 


1  1 


J_  J_ 
a" 
*^0 


exp  (- 


exp  (- 


1  (S-r,)2j 


1  (5-Rq2)^ 
2. 


(21) 


This  ratio  is  based  on  the  normality  assumption  of  the  estimated 
reliabilities.  Reducing  this  equation  and  taking  the  natural  logarithm 
yields  the  form  of  the  test  statistic  used  for  the  SPRTs  in  this  thesis: 


T  (R-Rq)  1  (R-R^)' 


2 

<y  r 


-  In  (op  )  +  In  (op  ) 

Rl  Kq 


(22) 


where  R  is  the  reliability  as  estimated  from  samples  using  MLE  of  k 
and  9,  Rg  is  the  reliability  under  the  null  hypothesis  and  R^  is  the 
reliability  under  the  alternate  hypothesis.  The  variances  of  R-j  and 
Rg  are  computed  using  the  Cramer-Rao  Lower  Bound. 

The  decision  boundaries  for  the  sequential  tests  are  developed 
as  follows: 


B  <  L(R^)/L(Rg)  <  A 

(23) 

In  8  <  ln(L(R.,)/L(Rg))  <  In  A 

(24) 

17 


and  since 


=  1n(L(R^)/L{RQ))  (25) 

In  B  <  <  In  A  (26) 

A  rainiraura  sample  size  of  10,  with  a  truncation  decision  occuring  at 
1.5  times  the  fixed  sample  size  is  employed  for  all  tests  in  this  thesis 
For  a  log-likelihood  test,  the  truncation  decision  is  to  reject  if 
Z^  >  0  and  to  accept  if  <  0  at  the  truncation  point.  Figure  3  is 

a  graph  of  a  typical  truncated  log  SPRT  when  the  minimum  sample  size  is 
10  and  the  test  is  truncated  at  1.5  times  the  fixed  sample  size  of  30. 
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Figure  3.  Truncated  Log  SPRT 
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Computerization 

The  first  step  in  the  computerization  process  was  to  develop  a 
program  to  estimate  the  shape  and  scale  parameters  for  the  two-parameter 
Weibull  distribution.  This  program  is  then  incorporated  into  the  main 
program,  to  perform  the  sequential  testing,  as  a  subroutine.  As 
mentioned  previously,  the  subroutine  to  estimate  the  parameters  uses 
the  bisection  method  to  find  a  root  to  equation  (15).  With  each  increase 
in  sample  size,  new  estimates  for  k  and  theta  must  be  calculated  using 
this  additional  information.  Subroutine  SOLVE  is  called  each  time  an 
additional  observation  is  taken  and  parameter  estimates  are  needed. 

From  within  SOLVE,  subroutine  EVAL  is  called  in  order  to  evaluate 
equation  (15),  since  this  equation  involves  summations  and  summations 
of  logs.  Values  of  THAT  (estimate  for  theta)  are  also  computed  in  EVAL 
for  each  estimate  of  k  (khat). 

Complete  listings  for  the  two  programs  (one  to  estimate  the  alpha 
error  and  one  to  estimate  the  beta  error)  are  given  in  Appendix  A. 

Also  illustrated  in  Appendix  A,  is  a  complete  variable  description  list. 

The  estimate  for  Alpha  (program  one)  is  performed  as  follows: 

1.  Input  a  sequence  of  Weibull  deviates  with  Hq  true  (line  470). 

2.  Run  1000  sequential  tests  (Monte  Carlo  sample  size). 

3.  Number  of  test  rejections/1000  is  the  estimate  for  alpha. 

4.  From  the  1000  test  runs,  an  average  number  of  items  tested  to 
an  accept  decision-Hg  true  (Aq),  and  an  average  number  of  items  tested 
to  a  reject  decision-H^  true  (Rq).  can  be  calculated. 
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5,  The  average  sample  size  is  a  weighted  sum  of  the  above  mentioned 
averages,  calculated  by: 

AVESAM  =  a  Rq  +  (l-a)AQ  (27) 

The  estimate  for  Beta  (program  two)  is  performed  as  follows: 

1.  Input  a  sequence  of  Weibull  deviates  with  true  (line  480). 

2.  Run  1000  sequential  tests. 

3.  Number  of  test  acceptances/1000  is  the  estimate  for  beta. 

4.  From  the  1000  test  runs,  an  average  number  of  items  tested  to 
an  accept  decision-H.|  true  (A-j),  and  an  average  number  of  items  tested 
to  a  reject  decision-H.j  true  (R-j),  can  be  calculated. 

5.  The  average  sample  size  is  computed  as  follows: 

AVESAM  =  (1-6)R^  +  6  A^  (28) 
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IV.  Results  and  Recommendations 


Results 

Appendix  B  contains  the  computer  results  of  the  various  Monte  Carlo 
simulations.  1000  Monte  Carlo  test  runs  are  used  for  each  case.  Computer 
Execution  time  varied  between  less  than  1000  to  over  6000  seconds,  with 
longer  times  being  associated  with  those  hypothesis  tests  in  which  the 
difference  in  Rq  and  R.|  was  less  pronounced,  and  the  output  risks  are 
smaller.  A  minimum  sample  size  of  10  is  used  for  each  test  with  the 
exception  of  indicated  departures  during  sensitivity  analysis. 

Table  B-1  shows  the  sequential  test  results  for  testing  the  null 
hypothesis  Hq:Rq=.9  against  the  three  alternatives  H.|  :R.|  =  .854,  H.|:R.|  =  .81, 
and  H.|  :R.|  =  .729.  Three  desired  (input)  alpha  and  beta  risk  level  pairs 
are  used  for  each  hypothesis  test,  (.1,  .1),  (.2,  .2),  and  (.1,  .2).  The 
criterion  used  to  judge  the  "goodness"  of  the  test  plans  is  how  close 
the  actual  alpha  and  beta  errors  are  in  relation  to  the  designed  alpha 
and  beta  output  risks.  Comparing  the  results  here  with  those  obtained 
by  Ballard  shows  noticeable  improvement,  especially  in  those  cases 
concerning  alternate  reliabilities  of  .81  and  .729. 

Sensitivity  analysis  was  conducted  to  determine  the  effect  of 
changes  in  the  minimum  sample  size  and  truncation  point.  One  test  case 
was  selected  in  which  the  minimum  sample  size  was  changed  from  10  to  15 
(Computer  resources  prevented  the  selection  of  additional  cases  to 
perform  sensitivity  analysis  on).  Table  B-2  indicates  that  this  change 
had  a  minimal  effect  on  the  test  plans.  There  is  only  minor  variation 
in  the  a  and  3  levels  accompanied  by  a  small  change  in  corresponding 


sample  sizes.  Also  shown  in  Table  B-2  is  a  change  in  the  truncation  point 
Here  the  minimum  sample  size  was  returned  to  10,  and  the  truncation  point 
was  moved  from  1.5  to  2.0  times  the  fixed  sample  test  size.  Again  it  is 
noted  that  this  change  had  only  a  minor  effect  on  the  test  plans. 

To  further  balance  the  a  and  g  levels,  the  boundaries  of  the  test 
region  can  be  shifted.  Results  in  Table  B-1  show  good  reason  to  shift 
since  the  alpha  errors  are  consistently  below  design,  while  the  beta 
errors  are  consistently  above  design.  Rather  than  shifting  the  boundaries 
it  is  much  easier  to  make  a  shift  in  the  test  statistic  to  allow  movement 
into  the  desired  test  region.  In  this  case  a  positive  shift  is  required 
in  order  to  raise  the  alpha  error  and  lower  the  beta  error  to  more 
closely  resemble  the  design,  and  thereby  accomplish  a  balance. 

The  test  of  hypothesis  Hq:Rq=.9  versus  the  alternative  H.|:R.|  =  .81 
was  selected  to  perform  these  shifts  on  (again  computer  resources  pre¬ 
vented  the  selection  of  additional  cases).  The  three  pairs  of  a  and  g 
risk  levels  were  tested  using  various  shifts  in  the  test  statistic. 

Table  B-3  shows  the  results  of  these  shifts.  These  results  are  consistent 
with  theory  in  that,  as  the  shift  increases,  the  alpha  level  increases, 
while  the  beta  error  decreases.  This  was  the  anticipated  and  desired 
result. 

Once  the  output  risks  have  been  balanced,  they  can  be  raised  or 
lowered  by  changing  the  design  (input)  risks.  In  this  case  the  output 
risks,  after  balancing,  were  above  the  desired  a  and  g  levels.  Table 
B-4  shows  the  results  of  a  shift  in  the  test  statistic,  along  with  input 
risk  levels  designed  to  bring  about  the  desired  results.  The  input  a 
and  g  levels  are  well  below  those  which  are  desired. 
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Recommendations 


Future  research  in  this  area  should  concern  the  development  of  a 
test  statistic  that  would  require  no  shifting  or  balancing.  Other 
methods  for  estimating  the  variance  for  the  estimated  reliabilities  may 
be  developed.  Estimates  for  the  parameters  k  and  e,  may  be  accomplished 
using  the  linear  regression  method  or  the  method  of  matching  moments, 
rather  than  the  Maximum  Likelihood  method.  Different  Alpha  and  Beta 
risk  levels  might  prove  to  be  a  worthwhile  area  of  investigation.  In 
addition  to  this,  different  tests  of  hypotheses  might  be  developed  in 
order  to  extend  these  test  plans. 
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Appendix  A 
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Program  one  -  Performs  sequential  tests  of  hypothesis  concerning  system 
reliabilities  to  estimate  the  Alpha  error. 

Program  two  -  Performs  sequential  tests  of  hypothesis  concerning  system 
reliabilities  to  estimate  the  Beta  error. 

Note:  The  IMSL  Library  must  be  attached 

Variable  Descriptions 


ACCEPT 

- 

number  of  test  accept  decisions 

ACCSAM 

- 

accept  decision  sample  size 

AHAT 

- 

estimate  for  Alpha 

ALPHA 

- 

input  Alpha  error  probability 

ANOT 

- 

average  number  of  items  tested  to  an  accept  decision 

AVESAM 

- 

average  sample  size 

BACK 

- 

last  section  of  equation  (15) 

BETA 

- 

input  Beta  error  probability 

BHAT 

- 

estimate  for  Beta 

DSEED 

- 

seed  for  random  number  generator 

FIXN 

- 

fixed  sample  size 

FXL 

- 

equation  (15)  evaluated  at  the  left  bound 

FXR 

- 

equation  (15)  evaluated  at  the  right  bound 

ITER 

- 

number  of  iterations  (bisection  method) 

K 

- 

actual  value  of  k  used  in  generation  of  Weibull  deviates 

KHAT 

- 

estimated  value  of  k 

L 

- 

Monte  Carlo  repetitions 

LB 

- 

lower  boundary  for  the  sequential  tests 
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LCOUNT  -  number  of  truncation  decisions 

LK  -  evaluation  of  equation  (15) 

MIN  -  minimum  sample  size 

N  -  sample  size  (current) 

NN  -  number  of  deviates  to  be  generated  in  GGEXN 

NTRUN  -  truncation  point  for  sequential  tests 

R(I)  -  exponential  random  deviate 

RO  -  reliability  under  the  null  hypothesis 

R1  -  reliability  under  the  alternate  hypothesis 

REJECT  -  number  of  test  reject  decisions 

REJSAM  -  reject  decision  sample  size 

RHAT  -  estimated  value  of  reliability 

RNOT  -  average  number  of  items  tested  to  a  reject  decision 

STDRO  -  standard  deviation  of  Rq 

STORl  -  standard  deviation  of  R^ 

SUMX  -  sum  of  the  Weibull  deviates 

SUMLX  -  sura  of  natural  logarithms  of  Weibull  deviates 

T  -  value  of  time  for  current  run 

THAT  -  estimate  for  theta 

TOL  -  tolerance  for  bisection  estimation 

UB  -  upper  boundary  for  the  sequential  tests 

A 

VARRO  -  variance  of  Rg 

A 

VARRl  -  variance  of  R-j 

X(I)  -  random  Weibull  deviates 

XL  -  left  bound  for  bisection  routine 

XM  -  mean  and  standard  deviation  for  use  in  random  number  generator 


29 


* 

-  ■ 


XR  -  right  bound  for  bisection  routine 

ZN  -  test  statistic 
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PL0.CH65000.T3000.  T810332,LUSS1ER,4462 
ATTACH, INSL.IDsLISRARy,SN«ASD. 

LIBRARY. IRSL. 

FTN5,ANSI»0. 

L60. 

lEOR 

C 

C  THIS  PRG6RAN  PERFORMS  A  SEQUENTIAL  TEST  TO  ESTIMATE 
C  THE  ALPHA  ERROR.  HO  IS  INPUT  AS  THE  TRUE  CONDITION 
C  OF  RELIABILITY,  AND  ALPHA  IS  ESTIMATED. 

U 

PR06RAM  SE8 
DOUBLE  PRECISION  OSEED 
DIMENSION  )((300),R(S) 

REAL  K,KHAT,L6 

DATA  RO, Rl.FITN, ALPHA, BETA, MIN, L/. 9,. 729, 14.,.!,. 2, 10,1000/ 

PfilNTt,  R0,R1,F1)(N,ALPHA,8ETA,MIN,L 

Ksl. 

T*U. 

DSEEDM4&2.D0 

NTRUN>l.S4FnN 

NN>1 

C 

C  UPPER  AND  LONER  BOUNDARIES  FOR  THE  SEQUENTIAL  TESTS 
C 

UB*L0S((1. -SETA) /ALPHA) 

LBsLfl6(BETA/(l. -ALPHA)) 

C 

LC0UNT»0 

ACCEPT=0. 

REJECTED . 

ACCSAM«0, 

REJSAH»0, 

C 

C  CALCULATION  OF  MEAN  AND  STANDARD  DEVIAION  TO  BE 
C  USED  AS  INPUT  TO  6SEXN 

C 

<M»-(TIIK)/LOfi(RO) 

C 

C  MONTE  CARLO  LOOP  OF  L  REPETITIONS 

C 

DO  200  1*1, L 
M«0 

40  CALL  S6EXN (OSEED, XM,NN,R) 

N»N+1 

JI(N)«R(1)UI1/K) 

C 

C  6ENERATE  MIN  NEIBULL  DEVIATES 

C 

IF(N.LT.HIN)  60  TO  40 
C 

CALL  SOLVE(X,KHAT,THAT.ITER.N) 

IF(ITER.6E.40)60  TO  40 
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ESTINATE  THE  RELIABILITY 


RHAT*E1(P(-((T/THAT)«IKHAT)) 

VARIANCE  AS  ESTINATED  BY  THE  CRANER-RAQ  LONER  BOUND 

VARRO«(ROU2)t(L0G(R0)tt2)t(1.109-.SI4t(LO8(-LO6(ROI))^.6O8t((LO6 

C(-L06(R0)l)tt2)l/N 

VARRl«(Rltl2)t(L06(Rl)tt2)t(1.109-.S14t(L08(-L08(Rl)))^608t(ILOe 

C(-L06(RI)IIU2))/N 

STDRO-VARROtl.S 

STORl-VARRltt.S 

CALCULATE  THE  TEST  STATISTIC  AND  PERFORM  THE  SEQUENTIAL  TESTS 

ZN* ( (RHAT-RO) <12) / (2IVARR0) - ( (RHAT-Rl ) «2) / (2<VARR1 ) -LOS ( STDRl ) ♦ 
CLOG(STORO) 

IF(N.SE.NTRUN)  60  TO  4S 
GO  TO  46 
LC0UNT*LC0UNT+1 
IF(ZN.6T.0.)60  TO  48 
IF(ZN.LT,0.)60  TO  47 
IF(ZN.LT.LB)60  TO  47 
1F(ZN.ST.UB)60  TO  48 
SO  TO  40 
ACCEPT»ACCEPT+1 
ACCSAN>ACCSAH*N 
60  TO  SO 
REJECTsREJECTM 
REJSAN^REJSAN^N 

INTERMEDIATE  RESULTS 

IF(I,E9.L/5)THEN 
AHAT»REJECT/I 
ANOT«ACCSAM/ACCEPT 
RNOT«REJSAH/REOECT 
AVESAM»AHATtRNOT+ ( 1 . -AHAT) tANOT 
PRINTS, ’HO:  R»  ’,R0 
PRINT!, 'HI:  R»  ’,R1 

PRINTS, ’INPUT  ALPHA,BETA  =  ’, ALPHA,’,’, BETA 
PRINTS, ’INPUT  HO  TRUE’ 

PRINTS, ’MONTE  CARLO  SAMPLE  SIZE  M 
PRINTS, ’NUMBER  REJECT*  ’, REJECT 
PRINTS, ’NUMBER  ACCEPT*  ’, ACCEPT 
PRINTS, ’TEST  ALPHA  ’,AHAT 

PRINTS, ’AVE.  NO.  ITEMS  TESTED  TO  AN  ACCEPT  DECISION*  ’,ANOT 

PRINTS. ’AVE.  NO.  ITEfIS  TESTED  TO  A  REJECT  DECISION*  ’,RNOT 

PRINTS, ’AVE.  SAMPLE  SIZE*  ’,AVESAM 

PRINTS. ’NO.  OF  TRUNCATION  DECISIONS  ’,LCOUNT 

ENDIF 

IF(I.EB.L/2)  THEN 


mM£j£cr/i 
AN0T»ACCSAI1/flCCEPT 
RNOT«REJSflN/REJECT 
AVESAK>AHAT tRNOT* 1 1 . -AHAT ) t ANOT 
PRINTt.’HO:  R*  ’,R0 
PRINTI.’HI:  R»  ’,R1 

PRINTI, 'INPUT  ALPHA, BETA  »  ’.ALPHA.’, ’.BETA 
PRINTt,’ INPUT  HO  TRUE’ 

PRINT!, ’NQNTE  CARLO  SAMPLE  SIZE  ’.I 
PRINTt, ’NUMBER  REJECT*  ’.REJECT 
PRINT!, ’NUMBER  ACCEPT*  ’.ACCEPT 
PRINT!, ’TEST  ALPHA  ’,AHAT 

PRINT!, ’AVE.  NO.  ITEMS  TESTED  TO  AN  ACCEPT  DECISION*  '.ANOT 

PRINT!. ’AVE.  NO.  ITEMS  TESTED  TO  A  REJECT  DECISION*  '.RNOT 

PRINT!, ’AVE.  SAMPLE  SIZE*  ’.AVESAM 

PRINT!. ’NO.  OF  TRUNCATION  DECISIONS  ’.LCOUNT 

ENDIF 

200  CONTINUE 

AHAT*R£JECT/L 
ANOT*ACCSAM/ACC£PT 
RNOT*REJSAM/REJECT 
AVESAM*AHAT!RNOT^  < 1 . -AHAT) !ANOT 
PRINT!, ’HOs  fi*  ’,R0 
PRINT!, ’Hi!  R»  ’.Rl 

PRINT!, ’INPUT  ALPHA, BETA  *  ’.ALPHA,’, ’.BETA 
PRINT!, ’INPUT  HO  TRUE’ 

PRINT!, ’MONTE  CARLO  SAMPLE  SIZE  ’,L 
PRINT!, ’NUMBER  REJECT*  ’.REJECT 
PRINT!, ’NUMBER  ACCEPT*  ’.ACCEPT 
PRINT!, ’TEST  ALPHA  ’,AHAT 

PRINT!, ’AVE.  NO.  ITEMS  TESTED  TO  AN  ACCEPT  DECISION*  ’,ANOT 

PRINT!, ’AVE.  NO.  ITEMS  TESTED  TO  A  REJECT  DECISION*  ’.RNOT 

PRINT!, ’AVE.  SAMPLE  SIZE*  ’.AVESAM 

PRINT!, ’NO.  OF  TRUNCATION  DECISIONS  ’.LCOUNT 

STOP 

END 

SUBROUT I NE  SOLVE  U , KHAT , THAT , I TER . N) 

C 

C  SUBROUTINE  TO  ESTIMATE  THE  SHAPE (KHAT)  AND  SCALE 
C  (THAT)  PARAMETERS  OF  THE  TNO-PARANETER  NEIBULL 
C  DISTRIBUTION,  USING  THE  BISECTION  METHOD.  EVAL 
C  IS  CALLED  IN  ORDER  TO  EVALUATE  THE  ESUATION  TO 

C  NHICH  A  ROOT  IS  REQUIRED. 

C 

REAL  KHAT 

TOL*. 00001 

ITER*0 

XL»0.1 

XR»5.0 

KHAT»XL 

CALL  EVALIX. KHAT, XL, XR.FXL.FXR, THAT, N) 

KHATrXR 

CALL  EVAL(X. KHAT, XL, XR.FXL.FXR, THAT, N) 

25  IF(FXLtFXR)  30,100,35 
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30  IFiJ(R-)(LLT.2.IT0L)  60  TO  100 

IFI1TER.6E.40)60  TO  100 

TENP'XL 
XL»()(L»)(R)/2 
KHAT«<L 

CALL  EVAL(X,KHAT,tL,l(R,FlL,FIR.THAT,N) 

ITER»ITER+1 
60  TO  25 
35  XR«XL 

XL*TENP 
<L»aL+«R)/2 
KHAT=XL 

CALL  EVAL(X,KHAT,XL,XR,FXL.FXRJHAT,N) 

ITER  *  ITER+1 
60  TO  25 

100  KHAT«(XL+XR)/2 

CALL  EVAL(X,KHAT,XL,XR,FXL,FXR,THAT,N) 

RETURN 

END 

C 

C  SUBROUTINE  TO  EVALUATE  EQUATION  (IS).  ROOT  TO 

C  THIS  EQUATION  6IVES  AN  ESTIMATE  FOR  K.  THIS 

C  EQUATION  REPRESENTS  THE  FIRST  PARTIAL  DERIVATIVE 
C  OF  THE  L06-LIKELIHOOD  FUNCTION  NITH  RESPECT  TO  K 

C  FOR  THE  TNO-PARAHETER  NIEBULL  DISTRIBUTION. 

C 

SUBROUTINE  EVALlX.KHATjXL.XR.FXL.FXS.THAT.m 

REAL  X(300),KHAT,LK 

SUNX>0. 

SUNLX»0. 

BACX>0. 

DO  SO  IM.N 

SUnX*SUNX>X(I)ttKHAT 

SUNLX<SUKLXK06(X(I)) 

THAT*(SUNX/I)tt(l/KHAT) 

BACK»BACK+  ( (X  ( I )  /THAT)  MKHAT)  ILOS  (X  ( I )  /THAT) 

LN» I /KHAT-I IL06 ( THAT) ♦SUHLX-BACK 
50  CONTINUE 

IF(KHAT.Eg.XL)FXL>LK 

IF(KHAT,EB.XR)FXR«LK 

RETURN 

END 
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PLl,Cn65000,T6000.  T8i0332,LUSSIER,4462 
ATTACH, IHSL, IO«L!BRARY,SN>ASO. 

LIBRARY, INSL. 

FTN5,ANS1»0. 

LSO. 

tEOR 

PR06RAH  SEQ 
C 

C  THIS  PROBRAH  PERFORNS  A  SEQUENTIAL  TEST  TO 

C  ESTIHATE  THE  BETA  ERROR.  HI  IS  INPUT  AS  THE 

C  TRUE  CONDITION  OF  RELIABILITY,  AND  BETA  IS 

C  ESTINATED, 

C 

DOUBLE  PRECISION  OSEEO 
DINENSION  l((300j,R(S) 

REAL  K,KHAT,LB 

REAOl , RO, R1 , FIXN, ALPHA, BETA,  HIN,  L 
PRINTI,  RO,Rl,FI)(H,ALPHA,BETA,HIN,L 
K»I. 

T»16. 

OSEED*4462.00 

NTRUNM.StFHN 

NN«1 

C 

C  UPPER  AND  LONER  BOUNDARIES  FOR  THE  SEQUENTIAL  TESTS 
C 

UB«LOS((l, -BETA) /ALPHA) 

LB»LOS(8ETA/(l. -ALPHA)) 

C 

LC0UHT»0 

ACCEPT'D. 

REJECT'D. 

ACCSAH'D. 

REJSAH'D. 

C 

C  CALCULATION  OF  HEAN  AND  STANDARD  DEVIATION  TO  BE 
C  USED  AS  INPUT  TO  BBEKN 

C 

)(H'-(TIW)/L06(R1) 

C 

C  NONTE  CARLO  LOOP  OF  L  REPETITIONS 

C 

DO  200  I'1,L 
N'O 

40  CALL  B6EIN(DSEED,1IN,NN,R) 

N«N+I 

X(N)'R(l)tt(l/IO 

C 

C  8ENERATE  NIN  NEIBULL  DEVIATES 

C 

IF(N.LT.NIN)  SO  TO  40 
CALL  SOLVEIX, KHAT. THAT, ITER, N) 

IF(ITER.SE.40)60  TO  40 


c 

C  ESTIttATE  THE  RELIABILITY 
C 

RHAT«EXP(-nT/THAT)UKHAT)) 

C 

C  VARIANCE  AS  ESTINATED  BY  THE  CRANER-RAO  LONER  BOUND 
C 

VARR0>(R0tl2) t (L06(R0) U2) 1(1. 109-.S14I(L06(-L0S(R0) ) )+.608t ( (LOB 
CI-L06(R0)))II2))/N 

VARRl*(Rltt2)t(L06(Rl)tt2)t(1.109-.314«(L06('L06(Rl}))^608t(IL06 

C(*L0e(Rl)))tt2))/N 

STDRO«VARROU.S 

STDR1«VARR1U.3 

C 

C  CALCULATE  THE  TEST  STATISTIC  AND  PERFORN  THE  SEOUENTIAL  TEST 
C 

ZN>((RHAT-R0)tt2)/(2tVARR0)-((RHAT-Rl)lt2)/(24VAfiRl)-L06(STDRl)^ 

CL06(STDR0) 

IFIN.SE.NTRUN)  60  TO  4S 
SO  TO  46 

45  LC0UNT>LC0UNT>1 
IF(2N, ST, 0.160  TO  48 
lF(ZN.LT.0.)6fl  TO  47 

46  1F(ZN.LT,LB)80  TO  47 

1F(ZN.8T,UB)80  TO  48 
SO  TO  40 

47  ACCEPT«ACCEPTM 
ACCSAN«ACCSAH*N 
SO  TO  SO 

48  REJECT*REJECTtl 
REJSAN>REJSAH+N 

C 

C  INTERNEDIATE  RESULTS 

C 

SO  IFd.EB.L/SITHEN 

BHAT-ACCEPT/I 
ANOT>ACCSAH/ACCEPT 
RNOT>REJSAN/REJECT 
AVESAH* ( 1 , *BNAT) IRNOT^BHATIANOT 
PRINT*, ’HO:  R»  ’,R0 
PRINT*, ’HI:  R«  ’,R1 

PR1NT*,’INPUT  ALPHA, BETA  «  ’.ALPHA,’,’, BETA 
PRINT*, ’INPUT  HI  TRUE’ 

PRINT*, ’NONTE  CARLO  SAHPLE  SIZE  ’,I 
PRINT*, ’NUHBER  REJECT*  ’.REJECT 
PRINT*, ’NUHBER  ACCEPT*  ’.ACCEPT 
PRINT*, ’TEST  BETA’.BHAT 

PRINT*, ’AVE.  NO.  ITENS  TESTED  TO  AN  ACCEPT  DECISION*  ’,ANOT 

PRINT*, ’AVE.  NO.  ITENS  TESTED  TO  A  REJECT  DECISION*  ’,RNOT 

PRINT*, ’AVE.  SAHPLE  SIZE*  ’.AVESAN 

PRINT*, ’NO.  OF  TRUNCATION  DECISIONS  ’.LCOUNT 

ENDIF 

lFlI.EB.L/2)  THEN 
BHAT*ACCEPT/I 
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AMT^ACCSAfl/ACCEPT 

RNOT*REJSAH/REJECT 

AVESAH* ( 1 . -BHAT I tRNOT^BHATIANOT 

PRINTI,‘HOi  R«  ’,RO 

PRINT!, ’Hli  R»  ',R1 

PRINTI, 'INPUT  ALPHA, BETA  «  ALPHA.’,’, BETA 
PRINT!, ’INPUT  HI  TRUE’ 

PRINT!. ’NONTE  CARLO  SAMPLE  SIZE  ’,! 

PRINT!, ’NUMBER  REJECT*  ’, REJECT 
PRINT!, ’NUMBER  ACCEPT*  ’, ACCEPT 
PRINT!, ’TEST  BETA’, BHAT 

PRINT!, ’AVE.  NO.  ITEMS  TESTES  TO  AN  ACCEPT  DECISION*  ’,AN0T 

PRINT!, ’AVE.  NO.  ITEMS  TESTED  TO  A  REJECT  DECISION*  ’,RNOT 

PRINT!, ’AVE.  SAMPLE  SIZE*  ’.AVESAM 

PRINT!, ’NO.  OF  TRUNCATION  DECISIONS  ’,LCOUNT 

ENDIF 

200  CONTINUE 
C 

C  FINAL  RESULTS 
C 

BHAT»ACCEPT/L 

ANOT*ACCSAM/ACCEPT 

RNOT*REJSAN/REJECT 

AVESAM* ( 1 . -BHAT) !RNOT*BHAT!ANOT 

PRINT!, ’HOj  R*  ’,R0 

PRINT!,’Hlj  R*  ',RI 

PRINT!, ’INPUT  ALPHA, BETA  *  ’, ALPHA,’,’, BETA 
PRINT!, ’INPUT  HI  TRUE’ 

PRINT!, ’MONTE  CARLO  SAMPLE  SIZE  ’,L 
PRINT!, ’NUMBER  REJECT*  ’, REJECT 
PRINT!, ’NUMBER  ACCEPT*  ’, ACCEPT 
PRINT!, ’TEST  BETA’, BHAT 

PRINT!, ’AVE.  NO.  ITEMS  TESTED  TO  AN  ACCEPT  DECISION*  ’,ANOT 

PRINT!, ’AVE.  NO.  ITEMS  TESTED  TO  A  REJECT  DECISION*  ’,RNOT 

PRINT!, ’AVE.  SAMPLE  SIZE*  ’.AVESAM 

PRINT!, ’NO.  OF  TRUNCATION  DECISIONS  ’,LCOUNT 

STOP 

END 

SUBROUTINE  SOLVE(X,KHAT,THAT,ITER,N) 

C  SUBROUTINE  TO  ESTIMATE  THE  SHAPE (KHAT)  AND  SCALE 
C  <THAT)  PARAMETERS  OF  THE  TNO-PARAMETER  HEI6ULL 
C  DISTRIBUTION,  USING  THE  BISECTION  METHOD.  EVAL 
C  IS  CALLED  IN  ORDER  TO  EVALUATE  THE  E8UATI0N  TO 
C  NHICH  A  ROOT  IS  REQUIRED. 

C 

REAL  KHAT 

T0L*.00001 

ITER*0 

XL*0.1 

XR*5.0 

KHAT*XL 

CALL  EVALIX.KHAT,XL,IR,FJL, FIR, THAT, N) 

KHAT*XR 

CALL  EVALa,KHAT,XL,IR,FXL,F]R,THAT,N) 
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25  IF(FXLIFXR)  30,100,35 
30  IF(]|R-l(L.LT.2.tT0L)  60  TO  100 

IF  (ITER.  BE.  40)  60  TO  100 
TENP'XL 
XL>(XL*XR)/2 
KHAT>XL 

CALL  EVAL(X,KNAT,XL,XR,FXL,FXR,THAT,N) 
ITER«ITER+1 
60  TO  25 
35  XR>XL 

XL*TEHP 
XL>(XL>XR)/2 
KHAT>XL 

CALL  EVAL(X,I(HAT,XL,XR,FXL,FXR.THAT.N) 

ITER  •  ITER+1 
60  TO  25 

100  KHATa(XL+XR)/2 

CALL  EVAL (X.KHAT, XL, XR, FXL, FXR, THAT,N) 

RETURN 

END 

SUBROUTINE  £VAL(X,KHAT,XL,XR,FXL,FXR.THAT,N) 

C 

C  SUBROUTINE  TO  EVALUATE  EQUATION  (IS).  ROOT  TO 

C  THIS  EQUATION  6IVES  AN  ESTIHTE  FOR  K.  THIS 

C  EQUATION  REPRESENTS  THE  FIRST  PARTIAL  DERIVATIVE 
C  OF  THE  L06-LIKELIN000  FUNCTION  NITH  RESPECT  TO 
C  K  FOR  THE  THO<PARANETER  NEIBULL  DISTRIBUTION. 

C 

REAL  X(300),KHAT,LK 
SUNX>0. 

SUHLX^O. 

BACK«0. 

DO  50  I>1,N 

SUnX«SUNXtX(I)ltKHAT 

SUNLX«SUNLXH06(X(I)) 

THAT»(SUHX/I)II(1/I(HAT) 

SACK»BACK*  ( ( X  ( I )  /  THAT )  UKHAT )  tL06  ( X  ( I  >  /THAT ) 
LK*I/KHAT-IIL08(THAT)+SUHLX-BACK 
SO  CONTINUE 

IF(KHAT.£B.XL)FXL*LX 

IF(KHAT.EQ.XR)FXR«LK 

RETURN 

END 
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Ag  -  Average  number  of  samples  to  accept  Hg  when  Hg  is  true 

Rg  -  Average  number  of  samples  to  reject  Hg  when  Hg  is  true 

Aj  -  Average  number  of  samples  to  accept  Hg  when  H^  is  true 

-  Average  number  of  samples  to  reject  Hg  when  H^  is  true 
AVG  Sample  Size  -  Weighted  average  of  Ag  and  Rg  or  A^  and  R^ 


equential  Test  Results 
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Point  =  1.5  times  fixed  sample  size 


Truncation  Point  =  tabled  value  times  fixed  sample  size 


Sequential  Test  Hesults  With  Shift 


Minimum  Sample  Size  =  10 

Truncation  Point  '1.5  times  fixed  sample  size 


Table  B-4 

equential  Test  Plans  With  Shift  an^ 


sample  size 
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Block  20.  (Continued) 

Rf^=.90,  is  tested  against  alternative  reliabilities  of  .35^. 
.310,  and  .729.  Three  pairs  of  alpha  and  beta  risk  levels  are 
considered  for  each  test  ((.l,.l),  (.2, .2),  and  (.1,.2).  A 
truncation  decision  for  the  sequential  test  is  made  at  1.5  times 
the  fixed  sample  test  size  for  the  same  conditions.  One  thousand 
?<!onte  Cairlo  repetitions  are  used  for  these  test  procedures. 


